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11.  abstract  (AMofnum  100 


Lascpinduccd  thermal  acoustics  (LITA)  is  a  promisinj:  optical  (oor  Rjvc  mmng  technique  for  pasdynamir  mea¬ 
surement.  The  \*'^*  nonlinear  process  is  a  sequen---  of  i»o  opio  acowstic  effocis,  rtecirostncnoo  aixt  *hsoq>tio«\'rapM< 
tnermaliration,  and  the  acousto-optic  effect.  The  esoluiioncf  the  laser-induced  acoustic  structures  temporalis  modulates 
and  thereby  the  LITA  signal.  Time  resolution  s>f  the  signal  pros  ides  the  sound  speed,  thermal  diffusivity.and  acoustic 
damping  rate,  along  with  information  about  atomic  or  molciular  energy  transfer  rates.  LITA  can  also  measure  spectra  of 
both  the  real  and  imaginary  gas  susceptibility.  The  phssics  of  LITA  is  discussed  and  the  dcnsation  is  sketched  of  a  simple 
analytical  expression  that  accurately  describes  both  the  magnitude  and  time  history  of  the  LITA  signal.  Early  expcnmenul 
results  are  presented.  Sound  speeds  accurate  to  0  ''T  ami  trartspiHt  propemes  accurate  to  hase  been  measured  in 
a  single-shot  wiihoul  calibration  More  realistic  rruvleling  »binjlJ  dramaiicalK  impro'C  transport-propens  measurement 
LITA  "spectra”  base  been  taken  of  weak  spectrat  lines  of  \  "i  c one emr.v. ions  lcs»  than  50  ppb  Signal  retlcciisilies  as 
high  as  10"^  have  been  estimjiol  New  applu  of  I  11  v  >n,  luding  seliVimeus.  are 'uggesied 
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I.ASr-R-INDUCRD  IHERMAL  ACOUSTICS  (LITA) 

FOUR  WAVE  MIXINC.  MF.ASUREMENT  OF  SOUND  SPEED,  THERMAL  DIFFUSIVtTY.  AND  VISCOSITV 

B.R.  Cummingf 

ORArn  ATT  AMtosAunTAi.  Laboratories 

('AI  IKIRVIA  iNSTTTVTF.Of  TEniNOI/)OV 

PaMdcna.  rA9||25 

Laser- inducni  thermal  acr^usiKs  •(  ITAi  it  a  prrtmitinn  optical  (nur-«t-ate  mitiiif  iectHiM|ue  for  fardynamic  meiRMirtiwcm.  The 
V '  nonlinear  prixTtt  it  a  vnjuerKe  of  tw  opto-acoutiir  effecu.  eleclfounction  and  alMOtptionAipid'lhefmaiimiaR.  and  die 
AtHitio  iiptii  cffcii  The  ctitiuiionitf  the  later-induced  acowtiic  ttniciaires  temporally  modnlMc*  v***  and  Iherctryihe  LfTA  ugnal 
Time  reviluiHin  o(  the  ti|tnal  prot  idet  the  sound  tpord.  thermal  diffuutily. and  acoustic  danpuif  rate,  alonf  with  mformatlon  ohoul 
auenu  le  moiecuLar  enerf;>  iranifa  rates  I  ITA  can  afut  measure  Sfweva  of  both  the  read  and  unaRMiary  gm  stncepuNItly.  The 
physKt  (tt  I  ITA  i«  divusvd  and  the  dentation  it  Uetched  of  a  umpfe  analyiieal  etprettaon  that  acruraariy  dcactihet  both  the 
rnaRniiudeandtinKhiitortofthel  ITA  ti^rtal  Farit  espmgaenulrrsulu  are  presented  Soamd  speeds  tawm  lo0.5*t'  and  tranipon 
prtipertiet  a,. urate  u>  Kate  Krm  mrasurod  m  a  unftr  shot  ntthoui  raiibraiatn  More  maliaic  wodclimt  limld drvnaucaily 
iinf»otcirant{t,et  po>fvrit  mraturemeni  I  ITA  t|>n.ua'  hate  horn uReii of  »eaii  ipretrai  Imesof  N()>  mcoaccntraliotH lets  than 
Vlppb  Sicnal  rrlWvtit  met  at  hifch  at  >•  '  hat e  hreti  ettimaied  Sr»  yplw  itinni  of  I  rr>  mi  IuiImi([  clnrimriry  arr  lurrtfiiT 

IntroductRUi 

A  ric«  eftkaJ  Ju^rsHiK  irthoi^o.*  f.c  unfir  tfae  mraturrmeM  of  icaidysutMc  |uraw>e<ers  has  hem  detelo|ird  at 
•  •AU’li  U:it  u\hrit«4iir  oimoJ  I  at<r  Induced  Thermal  AuttstKS  il  fTA).  uses  the  Utne-rctofted  Uf!Aif  fcnenacd  *u  f<Mr-**»e 
mttini;  i<<  >Tieaturr  iK*  w«u?tj  tpord  thermai  diffusitut  aatd  aoustu  danpMif  coeflkscfM  (and  there  In  rncosny)  of  a  pas  If  a  pat 
hchatet  a^ieJin^  to  a  lrv>«n  ,.alies.  rvaataei  of  tute.  the  temfiename  of  the  pas  may  hr  wfemd  fma  the  lOMd  speed  Tins 
i<xhnM|tK  It  timilar  ti>  the  I  ll'S  ‘laKr  induced  (dionon  scjoermpt  leelMM^ne  for  ihr  meaiurrmem  of  ttnmd  ipeeds  m  rrytuH*  * 
I'sind  compfcielt  urtcafihraied  unfle  shot  l.fTA.  the  tempmaurr  of  the  lahormory  atr  nat  aaanud  snttaa  03^.  hmtied  Ky 
uncerta.nit  <<(  the  imdocular  «>ei|chi  of  the  Stmsd  m  j«td  heam  peotnevy  The  thermal  (MTusnsiy  wd  acowAK  tfampmp  cacfftocM 
«ere  meavured  •nhin  Vo;  (4  puNoVd  taiues.  limited  hs  hmie  heam  star  and  thcrmiliqlion  ?me  dkets.  nhtch  nttr  tfnorod 
('Oi'iNjiion  lo  tmpio^r  ihr  «cw»t  H  irw  trsttpcri  prrprns  ineamtienuai  wnoTtes  thf  dnerinmAieui  of  thermaluaiion  rfs.  and 
the  mraturerTtmi  <»f  tvan  pn'ftlet 

The  prsenrtrs  of  the  I  ITA  fsprtimrnt  is  sttmtarsothaiof  e'AWAmdurcmann  StehetIUnMn  tfirrmptor  DFWM  (dcpmenac 
fttur  tkjtr  mmnpi  Httueter  the  nnrtlmrjrRt  prodminpihe  I  ITA  upnaf  o  a cmdmaiion of  ofUfMeonsiac and acctmap-oCAtcefTeess 
me  the  ttfierrm  suprrpnitirei  of  lyujMnm  mothantcal  sAjaes  used  « ihr  iRhe?  trOsa^acs  Aa  aadytital  espretston  for  the  Rhanlmc 
rrp!;::»l.  ;»;«  inc  i  IT  A  u^sor  has  hrvn  denied  assnuunp  laww  hydrodynamic  pai  hduntcr*  TVs  analysn  has 

thimn  ihsi  I  (T  A  tifruii  mat  iitaf  I  II  ■  later  mduted  dnrre  tceme >  re  DTWAf  stftuls  m  tMrnsaiy  «hcn ceuapanHr  lanes  arr  used 

!.>  djM  ;«t->  iptit  asftdcts.  fdiemunnw  hate  pnsdmod  I  ITA  stpnds  TV  hn(  is  ctettrtwantW*.  ta  nVeh  nudccnics  of  a 
(stlari/jMe  pat  minimi/e  ihnr  paemtal  enerps  hs  c  Vsarrmp  m  reftems  of  Vfh  cfticaf  fViinc  hcM  laamitiy*  WVn  ttit  ofltcal 
field  It  remmed  the  c  lusaered  pas  molnukrs  are  rtteaard  em— n  iwand  unset  ta  Ve  manarr  of  a  popped  halincn.  EVrtwtarteTtcm 
produces  a  I  IT  A  upnaf  prepnrtanul  to  the  si^aare  of  dv  rod  pan  of  ihr  pas  mnrfadalRy  dua  may  hr  nonrcaonaaa  or  rettmandy 
etshaned  The  second  <t«o  jcraruK  {thfncuacnon  h  dmaiaam  IcfVand  hy  rapid  nentaduBtse  decay  or  'dttrmiliaaaaaa*  of 
ahs>«ttnd  enrrpt  Rapid  tfsrrmdtrjainii  cnatn  ripd  espmiti'n  of  tfu  pas,  taduemp  siamd  s»om  TV  LiTA  akpaul  fVm  dai  clltcl 
II  (irtp>«iM«uia»Utr  Mpurr  <f  chr  unapmars  panef  mr  pas  nmefVNVs  TV  tuo phentuarna  prodacr  LIT A  lapnalt  nsih  ddfircaa 
iicrururrt  aiVnirtp  soninpif  the  nsntrdtaitmtodir  I  ITA  sapnd  fmn  each  cflcitt  TVn.  a  tpocirum  of  hodi  Vc  real  and  tatapwary 
(urti '  1  iV  pas  suscepiiVltit  mas  hr  mramrrd  asrap  I  fTA  TVsrrmcipRr  uai  dewnnanard  ttsmp  nead  vansatraN near  ttd  sndwn 
:»  linrt  i<r  tV  tdv'ranet  ati  iprv  tet  h  ‘  *  'le^*  iha»  V*  ppo>  . 


Ills  >mfratKtn«  .c .  la  n  ht>,«>s«aRm  tame  n  jAet  ml  drprud  spm  htdl  <«  ttsnuptunan  (vepmars  of  a  VpMy  sidhiaimd 
(a>  !h»i  ssuars  1  IT  A  •maio  amotf  laser  duacnctiattt  m-m  uHKh  mtta  cn  midocttljr  Mmr  ucaAes  and  drfdcar  nHdoraV 
^»4e%»v  tm-ejsiksm  1  IT  A  Hmrtiii  feesn  lariuip  esetirtj  tuar  mrrft  paemhinp  a  meanidV  narractara  daa  hampers  LIF  aad 
,<hft  epre to •«> . Sfc Sms^sfcs  !  J-  A  iswjiuimii  m  hrt*  •  ••Mwa*  ad*taaaie  •  Vm  m.dotalia  3<*saatrs  and  pdaeVnp  rmo 
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are  high,  as  in  many  flows  of  practical  gas-dynamic  and  aerothennochemical  interest.  For  example,  high  molecular  density  is  a 
requirement  lor  the  rapid  chemical  reactions  needed  in  advanced  propulsion  systems.  The  hydrodynamic  lifetime  of  the  LITA 
interaction  suggesu  new  applications  of  the  techni^,  among  which  arc  measurement  of  velocity  by  analogy  with  laser  Doppler 
velocimetry  and  measurement  of  vorticity.  for  which  there  is  no  conventional  analog. 

This  anicle  outlines  the  development  and  early  experiments  of  this  measurement  technique.  First,  the  physics  of  LITA 
mcasurementi  is  discussed.  A  brief  sketch  follows  of  the  derivation  of  the  lauilylical  expression  for  the  LITA  signal.  The 
espcrimcntal  setup  used  at  (iALCIT  is  then  presented  along  with  the  results  of  the  early  experiments.  Finally,  future  directions  are 
suggcsicd  f(w  LITA  application  and  conclusions  are  drawn 

HlYSKS  Of  UTA 

I ITA  »»  j  laser  induicd  d>oamH  grating  technique  similar  to  CARS  and  DFWM.  The  principal  difTcrence  between  these 
irvhniques  and  LITA  is  the  i nonlinear)  process  that  produces  the  signal.  LITA  employs  a  comMnalion  of  opio-acoustic  and 
aiouvto  riftK  effects  Light  from  a  powerful,  fast-pulsed  "dover”  laser  (e.g..  1  •  lOfiO  mJ.  I  •  10  ns)  generates  acoustic  waves  by 
AiwtUK  rflocts,  ekctnntrKtion  and  riptd  ihermaliration 

IK*  i.emer  rllnt  o  ^  suted  b>  the  tendency  of  polan/jMe  molerules  U)  minimire  their  potential  energy  by  clustering  in  regions 
'!  high  t^SKai  eifxtm  held  inicnsiis '  FlectrosuKlion  has  the  effect  of  lowering  the  pressure  in  regions  of  high  field.  When  the 
firSJ  o  ifm.-»rvl  there  o  a  pressure  imhaUmr  lhai  propagates  as  acofiltf  wjivet  This  effect  is  proportional  to  the  real  part  of  the 

ti'  'wS.  .•ptiNh!* 


I V  ijrtft  rll.^i  H  a  rmiluptr  step  prnsess  that  relies  upon  irolecular  collisions  First,  mokculcs  of  the  gas  absorb  energy 
Miv-i  I.V  <Vi«rr  laser  and  suer  ii  in  etc  tied  stales  when  the  driver  later  resonates  tvith  an  absorption  line  of  the  gas.  Next,  inelastic 
’iesulw  wiiltuaeit  cttnsrn  part  of  the  ctcited-sute  energy  to  molecular  kinetic  or  rotational  energy.  This  eonversion  may  require 
e.;<  a  L'w  ;oitick'«»  te  g  l(i<  iv  mans  le  g  .  hundreds),  depending  upon  ihe  composition  of  the  gas.  Finally,  clastic  molecular 
.  .silisseoetltiihbraaeihe  seiocitirs  d  ihe  mokculrs  resulting  in  a  nse  in  lempcTaiure.  This  equilibration  usually  requires  only  a  few 
ir  g  h«r  to  teni  The  lerm  Dirrmali/aiion*  covers  the  conversion  of  laser  energy  to  thermal  energy.  Thermalixalion  is 
rviSs^xoal  »r'  the  tmagmars  part  d  the  gas  suHeptibility 

Inivdrr  fit  the  ihrmtali/jiion  usprrtdsace  acoustic  waseiefficieniiy.  ihe  drop  in  gas  density  accompanying  the  rise  in  temperature 
'e.uu  kharpls  lay  that  of  oquihhratMo  The  density  should  be  "frozen  '  on  Ihe  lime  scale  of  Ihermaliration  (rn).  Density  changes 
n  a  Rmhum  traiel  efurac  ten  site  ally  at  the  uumd  speed  t' .)  Thermaliratkin  must  occur  on  a  lime  scale  comparable  lo  or  less  than 
t.he  lime  !•«  densiis  changes  to  trasel  across  a  rharaciensiK  LITA  length  scale  <.\f).  the  length  scale  over  which  Ihe  inientily  of  Ihe 
>11  ft  la  t  t  arses 


Kreaxise  is  mserscly  pnqvutNmal  hs  the  collisKsn  rale,  this  effoci  becomes  more  cfficicni  as  density  increases.  At  atmospheric 
dmsiisrs  ar»d  letnperaiurrs.  1  ITA  data  showed  that  fast  thcrmalirjiion  times  were  of  the  order  five  to  fifteen  nanoseconds  for  the 
’•< )  rscued  states  used  in  these  rtpaiments  lof  the  order  ten  collitioni'*)  This  compares  with  •\g/ct  of  thirty  to  onc-hundred 
ruooseximds  The  acousiK  fields  genenaed  by  the  daver  laser  evolve  in  lime  and  space.  They  are  damped  by  ihcrmal  diffusiviiy 
I  iKiHiiy  I  nder  Ihe  condiiions  studied  in  these  expenmenls.  this  damping  occurs  on  the  microsecond  lime  scale. 


I  igM  fnwn  a  long-pulse  or  TW  "siiurce''  laser  Katiers  into  the  LITA  signal  by  an  acouslo-oplic  effect.  Because  the  susceptibility 
a  medium  is  nearly  linear  wiih  densiiy  hu  acousiK  disturbances  *.  the  density  perturbation  field  of  the  acoustic  waves  generates  a 
.(vrrspimdtftg  suv.cpiiNliiy  perturbaiKin  field.  The  suscepfiNliiy  field  scatters  light  from  Ihe  source  laser,  forming  the  LITA  signal. 
AivjIssm  of  the  iimehfsiory  of  the  Kuiicrcd  signal  provides  infiinnaiion  about  Ihe  speed  and  decay  of  Ihe  opio-acouslic  disturbances 
aivJ  ihrrrbs  the  srwmd  speed  and  transport  pmpenies  of  the  gat 


Ath^uig  a  linir-wasc  misifig  gainirir)  for  LITA  as  shown  in  Figure  I  (wovides  both  spatial  resolution  and  coherent  signal 
(ohantcmeM  The  driver  laser  is  spin  into  issrt  beams  which  cross  at  a  small  angle  (0.3  •  5  degrees)  in  the  gas.  The  LITA  signal 
ti  genrrjtod  <«ds  tn  the  sample  volume,  where  ihex  beams  inicnecl.  providing  spatial  resolution.  Because  the  driver  beams  are 
.leterrm  with  rath  laher.  they  interfere  with  each  other  in  the  form  of  an  ellipsoidal  grating  at  ihcir  intersection.  The  grating 
•  virlmelh  •«,  IS  gisen  by 


•  ».  f.  It  s»ii-  .Witer  laser  »  jtelengih  jnd  **  is  ibe  bea^  •  angle 
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f  igure  I  ^  li'rv*ar«t  vjiu*nng  I. ITA  geometry  wa.v  uvcit  in  these  c«pcrimcniv.  The  ilnver  beams  write  acoustic  gratings 
into  the  gao  IIk'  MHiric  beam  irKidcnl  at  the  Bragg  ancle  scatters  oil  the  acoustic  structures  into  a  coherent 
signal  Uvim 


The  opto  acousiu  ellixis  gi“ns*raie  asiHisiu  siruclures  (modulation  depih  Odl  '‘tI  initially  in  the  shape  of  the  cilipsoitfal 
grating  Tlte  nature  of  thi*  ak.ousiK  strus  tures  generated  by  the  tws»  effes  ts  is  difk'rent  l-lectrosiriction  forms  two  traveling  acoustic 
IS  jscpackets  iphotMHisi  and  a  lers  \seak  stationary  isobaric  density  grating  i  ihermon  >  each  with  the  wavelength  of  the  electric  field 
grating  Thermal i rat i<*«i  generates  ihe  same  acinisiic  structures  but  the  •  ihcrmon"  i\  dominant  The  phonons  counterpropagaie  in 
the  direction  parallel  lo  the  grating  se»  tor  and  tkxay  by  acoustic  damping  The  ihermon  decays  by  heal  conduction.  Light  scattered 
Iriwn  the  ssxiixe  laser  bv  the  jiousik  grating  structures  adds  coherently  when  the  M>urce  laser  is  incident  at  the  Bragg  angle  of  the 
grating  As  shown  in  l-igurc  I.  the  Bragg -scattered  light  emerges  in  a  lobe  iw  laser-like  beam,  which  may  b<  cBicienily  detected 
with  large  f-numher  optics  This  Icaiore  is  attractive  where  optical  access  is  timiieil  or  where  luminosity  mandates  large  f-number 
detection 


Hach  aciHisiu  siriKlurc  scatters  light  Iron,  the  s«Hirce  laser  onto  the  detex  lor  II  there  is  noconvixtion  vrl.x*ily.  light  scattered  by 
Ihe  photMxi  moving  upw;ud  in  l  igiirc  I  has  a  iXippIcr  "rcd-shifi  "  Light  ssatiered  bv  its  coniugair  phonon  has  a  Doppler  “blue-shift.” 
Light  scattered  bv  the  thernMin  has  no  Doppler  shill  1he  beating  between  the  light  ss altered  by  these  structures  modulates  the  signal 
seen  by  the  square-law  detector  Three  frequency  components  are  present  in  the  signal  modulation  (neglecting  damping);  a  DC 
comporrent.  a  comporK-nt  at  the  •Brillouin"  frequency.  •',!  Ay. caused  by  phonon/ihermon  signal  mixing,  and  a  component  at  twice 
the  Brillouin  frequenev.  causs'd  bv  iiiiving  of  signals  generated  by  the  coniugaie  phomms  A  strong  phononAhcmKNi  frequency 
ctrmponeni  is  only  present  if  Ihe  Ihermali/ation  effect  is  present.  This  additional  frequenev  component  allows  distinction  of  the 
thermaltralion  aiKl  elecirosiriclion  effects  when  the  two  are  of  comparable  magnitude  The  real  .and  imaginary  parts  of  Ihe  gas 
susceptibility  may  thus  K’  measured  independently  using  LITA. 


Analysis  ol  LITA  .Signals 


This  section  outlines  the  derivation  ol  an  analytical  expression  for  the  amplitude  and  lime  history  of  Ihe  LITA  signal.  This 
expression  allows  Ihe  magnitude  of  I.ITA  signals  to  he  estimated  «i  inmn  for  experimental  lca.sibility  studies  and  the  actual  LITA 
signal  In  be  interpreted  quaniiiaiiv  civ .  w  ith  very  g«KKl  experimental  agreement.  The  analysis  treats  the  opin-acouslic  effects  oitd  the 
hydrodynamic  evolution  ol  the  acoustic  structures  using  the  lineari/cd  hydrodynamic  equations  ' .  The  perturbation  solution  for  the 
far  field  of  light  scattered  off  a  dielectric  disiurbaiTcc  accounts  for  the  acousio-optical  effect''.  A  simple  analytical  result  is  obuiined 
for  short -duration  Gaussian  driver  beams,  inliniic-cxicnt.  monochromatic  source  beam,  and  singlc-ratc-dominalcd  thermal i7.ation‘\ 

Opto-acousiic  effects  appetir  as  source  terms  in  the  closed  set  of  lincori/cd  hydrodynamic  continuity,  momentum,  and  energy 
equations.  Thermal  i/at  ion  behaves  like  a  temperature  source  in  the  energy  equation  while  elccirnstriction  forces  the  momentum 
equation.  In  terms  of  perturbations  in  density  (/»).  temperature.  (7  ),  these  equations  ore 
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where  perturbation  terms  are  indicated  by  the  subscript  ‘1’,  =  V  •  «i,  and  =  d-ldx^  +  d'^/dip’  +  d'^ldz^.  The 

parameter  7  is  the  raUo  of  specific  heats  of  the  medium  (Cp/c„).  a  is  the  thermal  expansion  coefficient  {dp /dT)^),  c,  is 

the  iscntropic  sound  speed,  Dy  is  the  longitudinal  kinematic  viscosity  {l/p  {%,  +  (4/8)773)),  and  Dt  is  the  thermal  diffusivity. 
The  boundary  and  initial  conditions  are 

pi(r,<  =  0)  =  0;  V>l(r,<  =  0)  =  0;  ri(r,<  =  0)  =  0,  (6) 

Pi(r,t)->0;  0i(r,t)-^O;  Ti(r,<)-^0,  as  |r|  —  00. 

Here,  liie  mean  flow  velocity  is  assumed  to  be  zero.  This  assumption  leads  to  a  loss  of  generality  because  of  the  lack  of  Galilean 
invariance  of  the  system.  The  general  case  of  arbitrary  mean  flow  velocity  will  be  explored  in  an  upcoming  extension  of  this  work. 


The  electrostriction  source  term  (i/>e(r,  /))  is  approximately 

^  (7) 

P 

fnr  "ir.iall  density  perturbations.  The  parameter  H,i  is  the  total  energy  of  the  driver  beams  during  a  pulse,  P(t)  is  the  normalized 
temporal  profile  of  the  driver  laser,  /( r )  is  the  normalized  driver  laser  field  intensity  in  the  sample  volume.  This  term  is  proportional 
to  the  real  part  of  the  gas  susceptibility  at  the  driver  frequency,  u>,/,  where  \  =  r  /f  n  -  1  • 


The  thermalization  term  7)/,  ( r,  t )  describes  the  conversion  of  encigy  from  the  driver  lavr  beams  to  heat.  While  this  process 
may  be  arbitrarily  complicated,  a  single-rate  rate  equation  for  the  excited  state  energy  (f/(  r.  / ))  with  diffusion  effects  approximates 
the  character  of  the  process: 

7;/(r./)  =  ^f’(r./).  (8) 

Pnr, 

where 

+  ‘infhf  -  (9) 

(It 

r(rJ  =  0)  -  0. 

where  U(r,t)  is  the  excited-state  energy  density  field,  '><7,  is  the  rate  of  loss  of  cxcilcd-statc  energy  to  thermal  energy,  “inth  i**® 
rate  of  excited-state  energy  loss  other  than  to  thermal  energy,  Du  is  the  diffusivity  through  the  gas  of  excited  molecules,  and 
IS  the  heat  capacity  of  the  gas  at  constant  volume.  If  thermalization  is  rapid  eiKMigh,  /  / ^  behaves  like  a  delta  function  rendering 
the  hydrodynamic  behavior  thermal ization-raie  insensitive.  If  a  single  thermalization  rate  mode!  docs  not  adequately  describe  the 
system,  the  linearity  of  the  equations  allows  solutions  using  different  rates  to  be  superimposed  to  simulate  a  more  complicated 
energy  conversion.  Thermalization  was  approximated  by  independent  fast  (c.g.,  10  ns)  and  slow  (e.g.,  100  ns)  processes  to  fit  the 
LITA  data  presented  later. 

Equations  3.4,5,  and  9  may  be  solved  by  applying  Fourier  and  Laplace  transforms'  The  resulting  set  of  linear  algebraic 
equations  has  the  solution  in  terms  of  the  Fourier  transform  and  Laplace  variables  q  and  •«  rcspcruvely: 
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where  perfect  gas  behavior  has  been  assumed.  The  quantity  3/(7, .-)  is  the  characteristic  equation  of  the  system  of  algebraic 
equations.  ,  ,  ,  .  »  1  ,  . 

M((i.  s)  =  s'  +  ( 1  Di<r  )•'*'  +  (^;V  +  <  w  fhu’  )•“  +  fhu’-  n  i ) 

The  roots  of  the  cubic  equation  A/(  .*<.  f/)  =  0,  provided  .\y  is  much  largr'^  than  the  mean-free  path  ol  the  gas,  arc  very  nearly 
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where  P  is  the  classical  acoustic  damping  coefficient,  ( 1  /2[{  y-l)DT+  Dy]).  and  r,  is  the  iscntropic  speed  of  sound. 

A  partial  fractions  expansion  of  equation  10  in  the  roots  of  s  simplifies  Laplace  inversion.  There  are  four  terms  in  this  expansion. 
TWo  terms  corresponding  to  the  real  roots  of  the  denominator  combine  to  describe  the  formation  and  decay  of  a  thermon.  The  terms 
resulting  from  the  two  imaginary  roots  describe  the  formation  of  conjugate  phonons.  The  acousto-optical  effect  is  modeled  using 
the  perturbation  solution^  for  the  far-field  of  light.  scattered  from  a  source  beam  by  a  small  disturbance  In  the  susceptibility  at 
the  source-beam  wavelength,  x^i  ( r,  f ;  wq). 

F  it  * 

£;,(R,/;q)  =  --^^ex|»i(k,  •  P.  -  u;o/)\i(q./:‘*'o).  (13) 

where  q  is  the  change  in  wave  vector  from  the  source  to  the  scattered  beam.  The  scattering  direction  and  strength  are  dicuitcd  by 
the  spatial  Fourier  transform  of  the  susceptibility  disturbance: 

/'■St 

\l{q,/;^o)  =  J  r/*roxp(iq- (M) 
The  susceptibility  of  a  gas  is  closely  proportional  to  the  gas  density  for  small  density  changes,  thus. 
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The  right-hand  side  of  this  eouation  is  known  in  terms  of  the  driver  laser  parameters.  Thus,  combining  equations  10. 13.  and  15.  the 
scattered  light  field  is  known  in  terms  of  the  imposed  lasers  and  parameters  of  the  gas  medium  including  the  sound  speed,  thermal 
diffusivity,  and  acoustic  damping  coefficient. 


The  scattered  signal  emerges  as  a  coherent  beam  with  the  divergence  angle  limited  by  diffraction  from  the  small  sample  volume. 
Integrating  the  square  of  the  scattered  electric  field  over  the  detector  area  yields  the  total  LITA  signal  intensity.  Signal  tuength 
may  be  expressed  as  a  reflectivity,  or  ratio  of  the  signal  beam  to  source  beam  intensity.  The  magnitude  of  the  LITA  reflectivity  is 

of  the  order  10“^  for  room  air  with  a  thermal  modulation  depth  of  1%  in  a  two-millimeter-long  sample  volume,  without  resonant 

cnharKcment  at  the  source-laser  wavelength.  Experiments  at  GALCIT  used  a  ooe-Wait  CW  Ar  loo  User  for  the  source  laser.  Prom 
reflectivities  of  this  magnitude,  a  milliwatt-range  He-Ne  source  laser  could  provide  usable  signal  levels. 

In  the  limit  of  rapid  thcrmali/ation,  short  driver-laser  duratkm.  and  flat  source-beam  intensity,  the  expression  for  the  temporal 
LITA  signal  becomes 

till  \\n  ^  ♦jpxpj  ( Ifi) 

where  t/g  =  ‘irr /.\y  is  the  driver  beam  interference  grating  vector,  and  .1.  //.  and  ♦  are  simple  algebraic  functkms  of  gas  and  laser 
beam  parameters.  Where  the  approximations  listed  above  arc  not  valid,  spatial  and  temporal  convolutions  of  laser-beam  profilex 
augment  the  analysis.  While  generally  these  convolutions  require  numerics,  all  of  the  parameters  which  augment  the  analysis  arc 
readily  measured  except  the  thermal ixation  rates.  These  rates  may  either  be  inferred  from  the  LITA  data,  estimated  from  published 
(lata,  or  ignored  provided  they  arc  much  faster  than  the  acoustic  iiKNion. 


The  philosophy  behind  the  development  of  LITA  has  been  to  obtain  the  most  accurate  and  robust  gas-dynamic  measuremenu 
with  the  least  experimental  and  analytical  complexity.  Common  and  relatively  inexpensive  tuneable  or  fixed-frequency  lasers 
may  be  used  for  LITA.  While  LITA  requires  wdc-bandwidth  signal  detection  and  recording  and  geometrical  phase-matching, 
these  experimental  complications  are  offset  by  the  intensity  of  the  LITA  signals  and  the  ease  in  signal  analysis.  Furthermore, 
LITA  measurement  of  the  sound  speed  providci  the  gas  temperature  If  the  caloric  equation  of  state  of  the  gas  is  known.  Other 
spectroscopic  techniques  usually  measure  temperature  by  scanning  a  rotaiional  band  of  a  well-studied  diatomic  molocule  with 
accessible  transitions,  of  which  there  is  only  a  handful.  The  populations  of  these  bands  must  be  inferred  from  the  signals,  not 
a  trivial  practice  in  some  diagnostics,  and  then  fit  to  a  Bolt/mann  distribution  to  obtain  temperature.  LITA  thermometry  may  be 
done  entirely  nonresonantly  <k  resonantly  using  ntiy  absorption  line,  provided  some  absorbed  energy  is  thero^ii^  on  the  LITA 
timc-scalc  or  faster.  This  freedom  in  ahmrpiion-linc  choice  may  significantly  simplify  experimenu  by  eliminating  the  need  for 
frequency  conversions  into  the  IJV  or  IR. 


Figure  2  shows  a  schematic  diagram  of  the  experimentaJ  apparatus  used  in  LITA  experiments  at  GALCIT.  The  driver  laser  is 
a  home-brew  narrowband  pulsed  dye  laser  oscillator/amplifier  with  a  noise  bandwidth  of  about  3  GHz.  Pumped  by  a  Continuum 
Q-switched.  frequency-doubled  NdrYAG  laser  model  YG-660  at  about  150  mJ/pulse,  the  dye  laser  output  is  about  30  mJ  in 
approximately  5  ns  pulses.  The  laser  operated  from  about  587  nm  to  about  592  nm  during  these  experiments. 


Figure  2.  The  LITA  experiments  conducted  at  GALCIT  use  a  Nd:  YAG  pumped  dye-laser  for  the  driver  laser  and  a  CW 
argon  ion  laser  for  the  source.  The  LITA  signal  is  delected  by  a  fast  PMT  and  recorded  by  a  fast  DSO. 


Approximately  ihrcc-millimcier-wide  driver  beams  arc  formed  from  the  dye  laser  emission  by  a  50%  beam  splitter  and  a 
mirror  ni  \ .  The  crossing  geometry  of  the  beams  ensures  that  both  beams  travel  similar  lengths  to  the  sample  volume  so  that  phase 
coherence  between  the  beams  is  not  lost.  Mirrors  nit  and  i)i,\  arc  adjusted  so  that  the  driver  beams  arc  parallel  and  enter  the  400 
mm  lens  symmetrically.  The  driver  beams  intersect  at  the  focus  of  the  lens  where  they  approximate  finite  plane  waves  alxsut  200  fim 
in  diameter.  The  beam  crossing  angle  is  set  by  the  distance  between  the  parallel  beams,  adjusted  using  micrometer  translation  stages 
under  mirrors  iiii  and  ni  BOXCARS  geometries'*  with  crossing  angles  ranging  from  1.0®  to  2.3®  were  used  in  these  experiments. 
E)ctcction  of  the  driver  beam  pulse  by  a  silicon  PIN  photodiode,  p]  (Thor  labs,  model  DET-2SI).  triggers  the  acquisition  of  the 
LITA  signal.  The  source  laser  used  in  these  experiments  is  a  CW  argon  ion  laser  operating  at  488  nm  with  a  I  W  output  power 
(Spectra  Physics  Model  165).  The  .source  laser  should  have  a  coherence  length  longer  than  the  interaction  region,  typically  two  to 
ten  millimeters,  because  pha.se  noise  of  the  source  laser  manifests  itself  in  the  washing-out  of  the  modulation  of  the  signal  beam. 
Future  experiments  arc  planned  using  a  (lashiainp-pumped  dye  laser  with  a  pulse  duration  of  about  800  ns.  a  pulse  energy  up  to  3 
J,  and  a  linewidth  of  about  0.(KI.3  nm.  Under  the  same  conditions,  this  laser  will  provide  one  million  limes  the  signal  obtained  from 
the  Ar  ion  la.scr.  providing  incohcreni-scaitering-linmcd  sensitivity. 

The  source  beam  passes  through  a  two-lens  system  which  adjusts  the  beam  width,  ranging  from  about  200  /im  to  2  mm  in 
these  experiments.  The  mirror  in  \  points  the  source  beam  at  the  sample  volume.  Phase-matching  adjustments  are  made  using  a 
micrometer-driven  translation  stage  under  ii»  | .  A  silicon  PIN  photodiode,  monitors  the  intensity  of  the  source  beam  during  the 
LITA  interaction. 

The  detection  system  consists  ol  a  last  photomultiplier  tube  (PMT,  Hammamaisu  model  OPTO-8)  with  at  least  a  5(X>  MHz 
signal  bandwidth  and  optical  filters  to  prevent  contamination  of  the  signal  by  otherwise  scattered  light,  luminosity,  room  light,  etc. 
The  filters  include  a  10  nm  bandpass  interference  filter  centered  at  488  nm  and  a  40  mm  lens/20  ;/m  pinhole  spatial  filter.  The 
LITA  signal  path-length  to  the  detector  is  about  2.5  m.  The  PMT  power  supply  (Fluke  4 1 2B)  was  carefully  filtered  to  eliminate  RF 
noise-sources. 

Data  are  recorded  using  a  1  CiS/s.  eight-bit.  two-channel  digital  storage  oscilloscope  (DSO)  (HP  545 lOA).  This  records 
typically  5(X)-io-2(XX)-sample  lime  histones  of  signals  from  the  PMT  and  p>.  mggered  by  the  signal  from  /»|.  When  needed, 
ensemble  averaging  eliminates  shot  noise  from  weak  Lll  A  signals.  Otherwise,  measurements  require  only  a  single  driver  laser 
pulse.  Data  from  the  D.SO  is  down  loaded  via  a  GPIB  interface  to  a  computer  for  analysis. 


•iPBir 


During  the  earliest  experiments,  LITA  signals  were  measured  from  the  laboratory  air,  which  had  essentially  the  same  composition 
as  the  outside  air.  According  to  the  Southern  Cal  ifomia  Air  Quality  Management  District,  this  contained  less  than  50  parts-per-billion 
of  the  pollutant  NOo  during  the  days  of  experimentation.  In  spite  of  the  low  concentration  and  the  relative  weakness  of  NO2 
transitions  near  589  nm,  LITA  successfully  detected  the  trace  species.  In  later  experiments,  NO2  was  seeded  into  the  air  near  the 
sample  volume  by  reacting  copper  and  concentrated  nitric  acid  in  a  beaker.  Future  tests  are  planned  to  explore  the  parameters  of 
LITA  quantitatively  in  a  high-pressure,  heated  bomb  with  optical  access. 


Experimental  Results  and  Discussion 


The  earliest  set  of  LITA  experiments  using  laboratory  air  was  completely  nonresonant,  the  signal  arising  completely  from 
nonresonant  clcctrostriction.  In  these  experiments,  the  signal  scattered  from  the  Ar  ion  laser  was  very  weak,  roughly  30  photons  per 
nanosecond  (reflectivity  of  ~  10~^  with  the  1  W  source  beam).  Ensemble  averaging  provided  usable  signals.  An  average  of  256 
shots  yields  the  nonresonant  data  plotted  as  the  dotted  curve  in  Figure  3.  The  frequency  of  the  modulation  is  about  26  MHz,  twice 
the  Brillouin  frequency,  yielding  a  sound  speed  accurate  within  the  uncertainty  of  the  knowledge  of  the  laboratory  sound  speed 
(~0.5%).  The  signal  recorded  near  a  line  of  N  0  >  appears  as  a  dashed  curve  in  Figure  3  for  comparison.  Resonant  enhancement  of 
the  real  part  of  the  susceptibility  of  the  gas  is  responsible  for  this  signal  increase.  The  signal  recorded  near  the  peak  of  the  absorption 
line  is  the  solid  curve  plotted  in  Figure  3.  On  this  curve  the  dominant  frequency  modulation  is  at  the  Brillouin  frequency,  caused  by 
interference  between  the  signals  scattered  by  the  phonons  and  the  thermon. 


Figure  3.  The  nonrcs('nani  electrosirictive  signal  (dotted  curve),  near-resonant  electrostrictive  signal  (dashed  curve) 
and  resonant  signal  (solid  curve)  are  plotted  for  comparison.  The  resonant  signal  modulation  is  mostly  at  the 
Brillouin  frequency.  Off-resonance  signal  modulation  is  mostly  at  twice  this  frequency. 


A  LITA  “spectrum''  of  the  susceptibility  of  the  gas  appears  in  Figure  4.  The  spectrum  was  taken  by  scanning  the  driver 
laser  between  test  runs  and  is  intended  to  suggest  the  capability  of  the  technique.  It  is  not  to  be  interpreted  quantitatively  because 
the  discrete  frequency  step  in  the  spectrum  was  about  12  GHz,  much  too  large  to  resolve  lines  of  the  gas.  A  spectrum  of  the 
complex  susceptibility  of  the  gas  may  be  obtained  by  analyzing  each  measurement.  Future  spectral  studies  are  planned  when  the 
laser-scanning,  data-acquisition/down-loading/data-processing  systems  arc  automated. 
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Figure  4.  This  LITA  “spectrum”  of  uace  NOo  consists  of  256-shot  ensemble  averages  of  LITA  signals  taken  approxi¬ 
mately  1 2  GHz  apart.  In  the  foreground  the  signal  comes  primarily  from  resonantly-enhanced  electrostriction; 
in  the  background,  the  signal  is  characteristic  of  resonance  (thermalization).  The  structure  in  between  consists 
of  unresolved  lines  because  of  the  large  frequency  step. 


Seeding  the  gas  with  NOo  provided  LITA  signals  that  were  several  orders  of  magnitude  suonger  than  the  signals  obtained  from 
room  air,  with  estimated  reflectivities  around  10~^.  True  single-shot  measurements  of  the  sound  speed  and  transport  properties 
were  taken.  A  sample  single-shot  LITA  signal  appears  in  Figure  5.  The  driver  beam  crossing  angle  was  1.04*^,  yielding  a  grating 
wavelength  of  32.3  /tm  and  a  Brillouin  frequency-ef  10.7  MHz.  At  high  concenuation  of  absorbing  species,  the  thermalization 
effect  swamps  the  elecuosuictive  effect.  Figure  6  shows  a  four-shot  average  of  the  LITA  signal  obtained  when  the  beam  crossing 
angle  was  widened  to  2.3  °  (grating  wavelength  of,  14.5  /im;  Brillouin  frequency  of  23.9  MHz)  Damping  of  the  waves  occurs  on 
a  time-scale  proportional  to  the  square  of  the  grating  wavelength.  In  Figure  6  the  theoretical  LITA  signal  is  plotted  with  a  dashed 
curve,  showing  excellent  agreement  with  experiment.  The  sound  speed  and  transport  coefficients  used  in  the  fit  are  published  values 
for  dry  air.  In  these  tests,  thermalization  occurs  on  a  time  scale  comparable  to  the  inverse  of  the  Brillouin  frequency.  Thus,  the 
simple  single-rate  model  of  thermalization  which  applies  when  such  rates  arc  poorly  resolved  is  replaced  by  a  model  consisting  of  the 
superposition  of  two  single-rate  thermalizations:  a  fast  <^10  ns  process  and  a  slow  100  ns  process.  The  two  rates  and  the  relative 
balance  between  the  processes  are  three  independent  fitting  parameters.  A  better  fit  to  the  data  could  be  obtained  by  adopting  a  more 
accurate  thermalization  model  and  sourcc-bcam-profilc  model  (a  top-hat  profile  is  assumed)  and  including  the  effect  of  source-beam 
phase  noise.  Future  experiments  will  seek  to  refine  these  models.  Because  of  the  resolved  unknown  thermalization  behavior,  the 
uncertainty  of  transport  property  measurements  is  high.  Experiments  conducted  at  high  pressures,  where  thermalization  rates  are 
much  faster,  arc  planned  to  confirm  the  accuracy  of  the  expression  for  the  LITA  signal  with  no  degrees  of  freedom  (in  the  limit  of 
rapid  thermalization  rate). 


Figure  5.  [a]  This  single-shot  LITA  signal  of  N  02-secded  air  was  taken  with  a  grating  wavelength  of  32.3  fm.  [b]  This 

four-shot  average  LITA  signal  of  se^cd  air  was  taken  with  a  grating  wavelength  of  14.5  /tm.  Damping-rates 
increase  as  the  inverse  square  of  the  grating  wavelength.  The  theoretical  fit  to  this  data  is  plotted  as  the  dashed 
curve,  showing  excellent  agreement. 


Future  Directions 

It  is  premature  to  gauge  at  this  early  stage  in  the  development  of  this  technique  what  its  future  capabilities  or  applications  might 
include.  However,  there  arc  several  interesting  applications  already  under  consideration.  The  measurement  of  transport  properties  of 
gases  is  very  important  in  several  fields  including  gas-phase  chemical  modeling  and  hypersonics.  To  the  author’s  knowledge,  LITA 
provides  the  most  direct  nonintrusive  measurement  of  these  properties.  The  measurement  of  sound  speed  provides  the  temperature  if 
the  gas  composition  or  caloric  equation  of  state  is  known.  While  this  is  not  a  direct  temperature  measurement,  the  considerable  ease 
with  which  the  measurement  is  performed  has  merit.  Furthermore,  the  “direct”  measurement  of  temperature  via  fitting  populations 
to  Boltzmann  distributions  often  requires  a  great  deal  of  spectroscopic  and  energy-transfer  knowledge.  LITA’s  ability  to  measure 
temperatures  accurately  from  the  very  first  signals  without  any  calibration  speaks  in  it’s  favor. 

While  the  measurement  of  these  gas-dynamic  parameters  has  been  discussed  throughout  this  paper,  there  are  several  other 
applications  that  are  evident.  LITA  has  provided  information  about  the  rates  of  thermal ization  of  excited-state  energy.  LITA  could 
in  the  future  become  another  window  in  the  study  of  molecular  and  atomic  energy  transfer.  The  ability  to  detect  and  measure  the 
complex  susceptibility  of  trace  species  suggests  that  LITA  could  be  employed  in  atmospheric  pollution  detection  and  spectroscopic 
studies. 

Another  application  of  LITA  is  llow-vclocity  measurement.  When  the  acoustic  structures  generated  using  LITA  convect  in  a 
flow,  the  light  scattered  off  them  generally  has  a  Doppler  shift.  This  Doppler  shift  may  be  measured  by  time-resolving  the  LITA 
signal  heterodyned  with  the  unscattered  source  laser.  The  strength  of  the  scattered  signal  in  LITA  makes  this  technique  practical. 
Preparations  for  future  experiments  in  LITA  vclocimetry  have  begun.  Yet  another  potential  application  of  LITA  is  the  measurement 
of  strong  vorticity.  The  phase-matched  scattering  angle  of  a  LITA  structure  immersed  in  a  fluid  with  local  rotation  will  similarly 
rotate.  Thus  the  direction  of  emergence  of  the  LITA  signal  could  be  used  to  infer  vorticity. 

In  all  likelihood,  the  pursuit  of  these  various  applications  will  point  toward  others.  In  the  near  term,  however,  the  analysis  of 
LITA  signals  and  techniques  of  LITA  measurement  will  be  refined.  Careful  parametric  studies  will  be  performed  to  validate  the 
analysis  and  find  the  limits  of  the  technique. 


Conciusions 


Lascr-induccd  thermal  acoustics  (LITA)  is  a  lime-resolved  four-wave  mixing  (FWM)  technique  for  gas-phase  measurement  of 
sound  speed  and  transport  properties.  It  is  different  from  conventional  tour-wave  mixing  techniques  in  that  the  optical  nonlinearity 
is  a  ‘■cquence  of  oplo-acouslic  and  acousto-optic  effects.  The  interaction  lifetime  is  dictated  by  the  damping  processes  of  acoustic 
attenuation  and  heat  transfer,  which  occur  generally  on  much  longer  lime  scales  than  dephasing  and  quenching.  Thus  the  LITA 
signal  may  be  resolved  in  lime  in  a  single  shot  using  a  fast  detector  and  fast  sampling  oscilloscope,  not  the  multiple-shot  optical 
delay  scanning  techniques  often  required  for  lime-resolved  FWM.  LITA  has  demonstrated  excellent  sensitivity  in  the  first  round  of 
experiments.  Detection  limits  of  weak  spectral  lines  of  NO2  in  atmospheric  air  of  less  than  50  ppb  were  observed  using  relatively 
weak  lasers.  LITA  sound-speed  measurements  were  accurate  within  the  uncertainty  of  the  calculated  laboratory  sound  speed 
(0.5%).  LITA  can  also  measure  the  complex  susceptibility  of  a  gas.  Frequency  scanning  experiments  enable  the  measurement  of 
susceptibility  spectra.  In  gas  seeded  with  N  O2,  reflectivities  of  the  order  10  ^  were  estimated.  In  both  magnitude  and  lime-history, 
the  LITA  signal  is  accurate' y  represented  by  an  analytical  expression  for  the  LI  FA  signal  derived  from  the  linearized  equations  of 
hydrodynamics  and  light  scattering. 
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